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Although the choice of a factorization scheme is as important as the choice of a factorization scale, 
the dependence of theoretical predictions (at finite order) on the choice of a factorization scheme 
has been little investigated. This is due to the fact that the freedom in the choice of a factorization 
scheme is enormous, even at NLO. One of the reason why to study factorization schemes is 



QJ the possible exploitation of the freedom in their choice in the construction of NLO Monte Carlo 

I event generators with NLO initial state parton showers. However, the ZERO factorization scheme, 

which should be optimal for such Monte Carlo event generators, has turned out to be practically 
^ inapplicable although it appears at first sight as reasonable. A detailed analysis has then shown 

that if some given NLO splitting functions do not satisfy a certain nontrivial condition, then the 
— . corresponding factorization scheme has some restrictions on its practical applicability. Relevant 

technical details of the discussed issues are the content of this short contribution. 
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1. The freedom in the choice of a factorization scheme at NLO 

To illustrate the issue of the freedom in the choice of a factorization scheme, consider a struc- 
ture function F(^x, Q^) which is expressed by the formula 

F{x,Q^)='£ f^Ci(-,Q\M,Fs)Diiy,M,FS), (1.1) 
f Jx y \y J 

where C;(x, 2^,M,FS) stands for the corresponding coefficient functions and D, (x,M,FS) repre- 
sents the parton distribution functions. The coefficient functions are fully calculable within the 
framework of perturbative QCD and can thus be expanded in powers of the QCD coupling param- 
eter a = Us/ 71 

Ci{x,Q\M,FS) =C;°^5(l-x) + a(At)C^^^(x,e^M,FS)+^(^^2(^)). (1.2) 
The parton distribution functions satisfy the evolution equations 

^^^-(«)E/'^«(^.M.FS)Z>.(,,..PS,. 

where the splitting functions Pij{x,M,FS) can be expanded in powers of a{M) 

Pij (x,M,FS) = Pjf^ (x) + a(M) P^p (x,FS) + ff{a^{M)) . (1.4) 

The higher order splitting functions pj^^ {x,FS), k> I, which can be chosen at will, can be used 
for labeling factorization schemes. At NLO, a factorization scheme is fully specified by the cor- 
responding NLO splitting functions P^j\x). If the series (1.2) and (1.4) ai^e summed to all orders, 
then the structure function f (x, Q^) is independent of the renormalization scale pL, of the factoriza- 
tion scale M and of the factorization scheme FS, however, finite order theoretical predictions for 
the structure function depend on these unphysical quantities. 

Changing the factorization scheme of the parton distribution functions is described by the 
formula 

A (x,M,FS) = A (x,M,FSo)+a(M)£ f ^ Tp^(-,FS,FS^ Dj{y,M,^^o) + ^{a^(M)) . 

j -Jx y \y / 

(1.5) 

The Mellin moments of the matrix function 7^^^' (x,FSi ,FS2) are determined by the following ma- 
trix equation 

T(i)(«,FSi,FS2),p(*'*(?i)j -ftT(i)(«,FSi,FS2) =P^^H",FSi)-P(')(?i,FS2). (1.6) 

The solution of the preceding equation can be expressed in analytic form, however, the Mellin 
inversion to x-space has to be calculated numerically in general. The transformation formula for 
the NLO coefficient functions reads 

4\x,Q',M,FS) =cf\x,e^Mo,FSo)+£cf (x)ln^ + 7)W(x,FSo,FS)) . (1.7) 
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2. The ZERO factorization scheme 

The so-called ZERO factorization scheme is defined by the condition that the NLO splitting 
functions Pj^j\x) vanish. The NLO initial state parton showers are thus formally identical to the LO 
ones in this scheme and therefore can be generated and attached to NLO QCD cross-sections by the 
standard algorithms. Using the ZERO factorization scheme, we can thus obtain NLO Monte Carlo 
event generators that involve NLO initial state parton showers without any change in the current 
algorithms. 

To obtain the coefficient functions and the parton distributions in the ZERO factorization 
scheme, we have to calculate the matrix function 7^.^^'(x,MS,ZER0) ^ This was performed for 
three and four effectively massless flavours. The obtained results are surprising because for x < 0. 1 

tI-I\x,MS,ZERO) ^ QjX-^ with (?if = 3) = 4.63 and («f = 4) = 3.85 . (2.1) 

The ZERO coefficient functions and the ZERO parton distributions diverge for low x in a. similar 
way as the functions 7^y''(x,MS,ZER0). The ZERO singlet parton distributions are plotted for 
X G (10^^, 10^') in Figure 1. The divergent terms cause problems in numerical calculations, and 
moreover, it is likely that the mutual cancellation of the divergent terms in expressions for phys- 
ical quantities, such as (1.1), is incomplete at NLO, which can lead to unreasonable theoretical 
predictions. From the practical point of view, the ZERO factorization scheme is thus in general in- 
applicable, but it can be proven that there are no problems with its applicability in the non-singlet 
case. 

3. The condition of applicability 

A given factorization scheme obviously has no restrictions on its practical applicability if the 
conxsponding coefficient functions, cross-sections at paiton level and parton distribution functions 
behave for low x in a similar way as the MS ones. If we consider only such factorization schemes 
for which the corresponding matrix P^^^ (x) has the same structure as the MS one, that is 

= W = 5uPi'Ax)+Pil^\x), P^lix) = P^^lix) = P^'Jix), 

(3.1) 

then a detailed analysis of the solution of the equation (1.6) shows that the preceding condition 
of applicability is fulfilled if and only if the Mellin moments of the corresponding NLO splitting 
functions Pil^^{n), P^i'^(n), Piy\n), P^cin), Pg^in) and P^^i{n) satisfy the equation 

P^in) {pS\n)-PSi{n)-b) [p^^J^n) -PiS{n,MS)) + 

+ {pS\n)-PSi{n) +b) -p(',)(n,MS)) - 

-2P(2(.)P£)(.) {Pii>'(n)+p!^(n) + 7.,pj^'(n) -F^in) - 

-P^P%,m) -P^(i%,MS) -2«fP^^)^(«,MS) +P^2(«'MS)) = (3.2) 

'According to the equation (1.6), the matrix function /^^'(x, ZERO, MS), which is required for the transformation 
of the parton distributions, is equal to — 7^.y(x,MS,ZERO). 
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Figure 1: Comparison of the ZERO and MS singlet parton distributions at M = 50 GeV. The MS distri- 
butions were obtained by evolving the starting distributions of the MRST98 set [1] with the fixed number 
of active flavours «f = 3 (only light flavours are taken into account). The ZERO distributions were calcu- 
lated numerically from the MS ones. Note that the gluon distributions are plotted in their absolute value 
because the ZERO gluon distribution is negative in the displayed region. The point where the ZERO gluon 
distribution changes the sign is close to x = 0.1. 

for n G ^ where 

^^^3 = {1.7329,4.6306}, ^,=4 = {1-7995,3.8458} , ^,,=5 = {1-9001,3.1798} . (3.3) 

4. Summary and conclusion 

It has been shown that not all NLO splitting functions that appear at first sight as reasonable 
specify practically applicable factorization schemes. The practical applicability of a factorization 
scheme is assured if the corresponding NLO splitting functions satisfy some nontrivial condition, 
which can be easily formulated in the space of Mellin moments. This condition is unfortunately 
not satisfied in the ZERO factorization scheme, which would otherwise be optimal for NLO Monte 
Carlo event generators. Hence, searching for a suitable factorization scheme which is close to the 
ZERO factorization scheme and satisfies the condition of applicability has already been started. 
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